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Abstract To find the discrete symmetries of a Hamilton operator H is of central 
importance in quantum theory. Here we describe and implement a brute force 
method to determine the discrete symmetries given by permutation matrices for 
Hamilton operators acting in a finite-dimensional Hilbert space. Spin and Fermi 
systems are considered as examples. A computer algebra implementation in Sym- 
bolicCH — h is provided. 
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1 Introduction 

In quantum mechanics the system is described by a self-adjoint (Hamilton) opera- 
tor H acting in a Hilbert space H. Here we consider the finite dimensional Hilbert 
space C n where H is a hermitian matrix [H El El HI [5] . One of the main tasks is to 
find the n x n unitary matrices U such that U*HU = H, where U* = C/ -1 . The 
n x n unitary matrices form the compact Lie group U(n). Note that if U*HU = H 
and V*HV = H then (UV)*H(UV) = H, where (UV)* = V*U*. Thus the set of 
matrices that keep H invariant form a group themselves [6]. 

An important finite subgroup of the group U(n) are the nxn permutation matrices. 
The number ofnxn permutation matrices is n\. For a given hermitian nxn matrix 

H we want to find all the nxn permutation matrices P such that 

P T HP = H 

where P" 1 = P T . These permutation matrices form a subgroup of all n x n per- 
mutation matrices. Obviously the nxn identity matrix I n satisfies I n HI n = H. 



Here we describe and implement in SymbolicC++ a brute force method to find 
these permutation matrices. 

After this finite group has been found one determines the conjugacy classes. Now 
for a finite group G the number of conjugacy classes is equivalent to the number of 
non-equivalent irreducible matrix representations. From the conjugacy classes and 
the permutation matrices we can construct projection matrices to decompose the 
Hilbert space into invariant sub Hilbert spaces |Bj. For example, if the permutation 
matrix P satisfies P 2 = I n , then n x = (J„ + P)/2, Tl 2 = (I n — P)/2 are projection 
matrices (with riill^ = 0) which can be utilized to decompose the Hilbert space 
C n into invariant subspaces. 

2 Examples 

We consider four examples: two Fermi systems and two spin systems. 

Example 1. Let c], Cj (j = 1,2,3) be Fermi creation and annihilation operators, 
respectively. Consider the Hamilton operator 

H = t(c\c 2 + 4 C 1 + 4 C 3 + 4 C 2 + 4 C 3 + 4 C + ^l4 C l + ^4 C 2 + hc\c 3 

and the number operator N = 4 c i+4 c 2+4 c 3- Then [H, N] = 0. Since [H, N) = 

we find N is a constant of motion, i.e. the total number of Fermi particles remains 

constant in the sense that if \n) is an eigenstate of the number operator N with 

eigenvalue n at time 0, then \n)(t) = e~ tHt l h \n) remains an eigenstate of iV with 
eigenvalue n for all times. Given a basis with two Fermi particles 

44lo>, 4410), 44|o). 

Then we find the matrix representation of H 

I k x + k 2 t -t \ 

H = I t h + k 3 t 

\ -t t k 2 + k 3 J 

The matrix representation of iV is the diagonal matrix 21%, where J3 is the 3x3 
identity matrix. For k\ 7^ k 2 , k\ 7^ k$, k 2 7^ k% no non-trivial symmetry is found. 
Also for k\ 7^ k 2 , k\ 7^ k 3 , k 2 = k 3 no non-trivial symmetry is found. For k = k\ = 
k 2 = k 3 we obtain the permutation matrix 

/0 1\ 
P= 1 . 

\i 0) 



Thus using the projection matrices 111 = (J3 — P)/2, n 2 = (I3 — P)/2 the Hilbert 
space C 3 can be decomposed into invariant subspaces. For the case k — k± — h% — 
ks we find the eigenvalues 2k — 2t and 2k + t (twice). 

Example 2. Consider the Hamilton operator (two-point Hubbard model) 

H = t(c{ t c 2t + c ] n c 2 i + 4 t e lt + c ] n c n ) + U {n n n n + n 2t n 2 J 

where := c^c^, riji := c^c,^. The operators c^, c]^, c^, are Fermi opera- 
tors. The Hubbard Hamilton operator commutes with the total number operator 
N and the total spin operator S z where 



2 1 2 

N := 2^(c] t c it + c] lCjl ), S z := - ^( c jt c it - cj 



We consider the subspace with two particles N = 2 and total spin S* z = 0. A basis 
in this four dimensional Hilbert space is given by 

C lt C ljJ0)' C lt C 2l|0)' C 2t C ll|0)' C 2t C 2i|0) ' 



We find the matrix representation of H for this basis. Using the Fermi anti- 
commutation relations and Cj||0) =0, CjjjO) =0 for j = 1,2 we obtain the matrix 

representation of H with the given basis 



H 



(U t t 0\ 

toot 
toot 

\0 t t U ) 



We find the four permutation matrices Pq — I4 = J 2 * I 2 — I2 <S> I2, 



1 

1 



P-2 



1 

1 



*/ 2 , Pz 
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0/ 


1*0 
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where we define the star product * of the two 2x2 matrices A, B as [B] 



«n ai2 \ ( b n b 12 

Q21 Q-22 / \ &21 ^22 



/an a 12 \ 

b n 6 12 

&21 &22 

\ a 2 i a 22 / 



We note that the star product of two 2x2 permutation matrices is a 4 x 4 permu- 
tation matrix. Here Pi is the swap gate. The four permutation matrices Pq, P\, 



P 2 , P 3 form a commutative group with P 2 = I4 for j = 0, 1, 2, 3. If P is an n x n 
permutation matrix with P 2 = I n then 



ni = -(/ n + p), 



n 2 = -(4-P) 



are projection matrices with nill2 = n , where n is the n x n zero matrix. Using 
the permutation matrices Po an d P3 (which form a subgroup) these projection 
operators can now be used to find the invariant subspaces 



/1\ 





V 1 / 



1 
1 

W 



/ 1 \ 





V- 1 / 



/ \ 

1 

-1 
V / 



These four vectors (after normalization) are known in quantum computing as the 
Bell basis [TJ [2J El HI [2] . This leads to the two invariant sub Hilbert spaces 



^(4/ lt |o) + 4 ; 4 t |o)), ^(4^10) + 4/ lt |o)) 



- 4Ai°»' 7f ( c l4 t l°) - 4At\°)) 

Example 3. Let a±, 02, (J3 be the Pauli spin matrices 



0"! 



1 

1 



a 2 :-- 



-i 

1 



0-3 



1 
-1 



Consider the Hamilton operators j7] 



H = huJias <g> I 2 + fiw 2 h ® cri + e(a 3 ® Gi) 



K = fkuia 3 (2) J 2 + fko 2 l2 + e(cri ® (T3) 

where for the second Hamilton operator K the interaction term is swapped around, 
i.e. CT3 <S> o"i — > 0\ ® cr 3 . This provides symmetry breaking. For the Hamilton 

operator H we find the symmetries (permutation matrices) 

Po = h, Pi = / 2 ©ai, P 2 = ai©/ 2 , P 3 = (Ti©(Ti 

where © denotes the direct sum. The four matrices form a commutative group 
under matrix multiplication. All satisfy P 2 = I4. Thus we can use the projection 



matrices 111 = (J 4 + Pj)/2, il 2 = (/ 4 — Pj)/2 to decompose the Hilbert space into 

two invariant subspaces. On the other hand for the Hamilton operator K we only 
find the identity matrix Pq = I4, i.e. no non-trivial symmetry is admitted. 

Example 4. Consider the Hamilton operator for triple spin interaction 

H = o\ © 02 © 03 . 

The eigenvalues of this hermitian and unitary 8x8 matrix are +1 (four-fold de- 
generate) and —1 (four-fold degenerate). Owing to this degeneracy one expects a 
"large" number of symmetries. Applying the SymbolicC+-|- code we find 24 permu- 
tation matrices listed Pq, Pi, ... , P23 with P = I 8 . They form a non- commutative 
group under matrix multiplication and are a subgroup of the group of 8 x 8 per- 
mutation matrices. We note that the Kronecker product © and the direct sum © 
of two permutation matrices is again a permutation matrix [8]. Now we can list 
the ones with P? = J 8 

Pi - 

P2 - 

P 5 - 

Pe - 
Ps - 

Pvt - 

Pis - 
P23 ~- 

The other ones can be found by multiplication of these permutation matrices, for 
example P3 = P1P2 etc. Thus the 24 matrices form a subgroup of the permutation 
group of 8 x 8 matrices. 

Another spin Hamilton operator studied is [9] 

4 4 

H = aJ2 *s0>30' + 1) + bj^ ^i(i) 
i=i i=i 

with cyclic boundary conditions, i.e. (7^(5) = cr 3 (l). Here a, b are real constants 
and 0i, a 2 an d 03 are the Pauli matrices. Thus the underlying Hilbert space is 



We have 



1 

1 

(h*(7i) © (01 *h 
1 





(cri *I 2 ) © (h*cri 
h ® o\ ® <Ti 


x° 1 
01 © J2 <8> ^2 
01 © 01 © 01 . 



1 

1 



01 




1 



1 




C 16 . Recall that 

a k {l) = a k ® h® h® h, cr&(2) = h ® Ok ® h ® h 

a k (3) = I 2 ® h ® cr k ® h-, c & (4) = h ® h ® h ® o k 
where k — 1, 2, 3. We obtain the symmetric 16 x 16 matrix for H 



/4a 


b 


b 





b 











b 




















o \ 


b 








b 





b 











b 




















b 








b 








b 











b 




















b 


b 














b 











b 














b 














b 


b 

















b 














b 








b 


-4a 





b 

















b 














b 





b 








b 




















b 














b 





b 


b 


























b 


b 


























b 


b 





b 














b 




















b 








b 





b 














b 

















b 





-4a 


b 








b 














b 

















b 


b 














b 














b 











b 














b 


b 




















b 











b 








b 








b 




















b 











b 





b 








b 


V o 




















b 











b 





b 


b 


4a / 



The Hamilton operator H admits the C± v symmetry group. The order of this 
non- commut at ive group is 8. One finds the following set of eight symmetries [H] 



E : 


(1,2,3,4) (1,2,3,4) 


C 2 : 


(1,2,3,4) -> (3,4,1,2) 


c 4 


(1,2,3,4)^(2,3,4,1) 


c\ 


: (1,2, 3, 4) -+(4, 1,2, 3) 




(1,2,3,4)^(2,1,4,3) 


i 

a v : 


(1,2,3,4)^(4,3,2,1) 




(1,2,3,4) (1,4,3,2) 




(1,2,3,4)^(3,2,1,4) 



which form a group isomorphic to C± v . The symmetries can be found by calculating 
the 16 x 16 permutation matrices such that H = P T HP. 

3 Code Description 

Algorithms for finding all permutations of a sequence of objects are described by 
Knuth [10] . For a given n the permutation matrices are generated with the follow- 
ing algorithm. The algorithm implements the nested loops 



For jo = 0, 1, . . . , n — 1 do 
For ji = 0, 1, . . . , n — 1 do 



For j n _i = 0, 1, . . . ,n - 1 do 
If Jo ^ ji ^ ■ ■ ■ ^ Jn-i then 

use the permutation (0, 1, 2, . . . , n - 1) ->• (j ,ji,j 2 , ■ ■ • , jn-i)- 
End loop 

End loop 
End loop 

Algorithm to find all permutation matrices. 

1. Create an array (jo,ji, ■ ■ ■ ,j n -i) of loop variables. 

2. Initialize jk '■= —1 for k — 0, 1, . . . , n — 1. 

3. Initialize the loop variable index i to i := 0. 

4. While z > 

(a) Iterate. 

Set ji := ji + 1. 

(b) Termination condition. 

If ji = n terminate this loop: 

i. Set ji := — 1. 
ii. Exit the nested loop. 



iii. GotoH 

(c) If jk = ji for some k £ {0, 1, . . . , i — 1} then goto HI 

(d) Enter the next nested loop. 
Set i :— i + 1. 

(e) Innermost loop completes. 

If i = n then use the permutation (0, 1, 2, ... , n— 1) — > (jo, ji, j2, • • • , jn-i) 
i.e. the permutation matrix P is given by 



Set i :— i — 1. 




The SymbolicC++ program [TT] utilizes the vector class of the Standard Template 
Library. The Hamilton operator refers to example 2 in the text (two point Hubbard 
model). 

// permutation. cpp 

#include <iostream> 
#include <vector> 
#include "symbolicc++ .h" 
using namespace std; 

int total; 
Symbolic H; 

void commutes (const Symbolic &P) 
{ 

if(P*H==H*P) cout « "PC" « total++ « "] = " « P « endl; 

} 

void f ind_perm(int n,void (*use) (const Symbolicfe)) 
{ 

int i , k ; 
Symbolic P; 
vector<int> j(n,-l); 

i = 0; j [0] = -1; 
while (i >= 0) 
{ 

if (++j [i] ==n) { j [i — ] = -1; continue; } 
if(i < 0) break; 

for(k=0;k<i;++k) if ( j [k] ==j [i] ) break; 

if(k!=i) continue; 

++i; 

if (i==n) 
{ 

P = 0*Symbolic("" ,n,n) ; 
for(k=0;k<n;k++) P(k,j[k]) = 1; 
use(P) ; 
— i; 
> 

} 

} 

int main (void) 
{ 

using SymbolicConstant : : i ; 
Symbolic sqrt2 = sqrt (Symbolic (2) ) ; 



Symbolic U("U"); Symbolic t("t"); 

H = ((U,t,t,Symbolic(0)) , (t, Symbolic (0) ,Symbolic(0) ,t) , 

(t, Symbolic (0) ,Symbolic(0) ,t) , (Symbolic(O) ,t,t,U)) ; 
cout « "H = " « H « endl; 
total = 0; 

f ind_perm(H.rows() .commutes) ; 
return 0; 

} 

A Maxima implementation is available from the authors. 

4 Conclusion 

We applied a brute force method to find all possible permutation matrices that 
provide symmetries for given Hamilton operators in a finite dimensional Hilbert 
space. With growing size of the Hamilton operators matrix representation finding 
the permutation matrices becomes very time-consuming. A more efficient approach 
would be to find only the generators of the group of permutation matrices that 
provide symmetries for a given Hamilton operator. Another open question is how 
this method can be extended to find other classes of symmetries. 
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